Resonant finite-size impurities in graphene, unitary limit and Priedel oscillations 
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Unitary limit for model point scatterers in graphene is known to reveal low-energy resonances. The 
same limit could be achieved from hybridization of band electrons with the localized impurity level 
positioned in the vicinity of the Fermi level. The finite size defects represent an easier realization of 
the effective unitary limit, occurring when the Fermi wavelength induced by the potential becomes of 
the order of the size of the defect. We calculate the induced electron density and find two signatures 
of a strong impurity, independent of its specific realization. The dependence of the impurity-induced 
electron density on the distance changes near resonances from oc r~ 3 to oc r~ 2 . The total number 
of induced particles at the resonance is equal to one per degree of spin and valley degeneracy. The 
effects of doping on the induced density are found. 

PACS numbers: 73.22.Pr,73.20.At,73.22.Dj 



I. INTRODUCTION 

In conventional three-dimensional electron systems 
probe charges are screened exponentially with the 
distance^. In degenerate metals at distances exceeding 
the screening radius the non-monotonic power-law tail 
develops, oc cos(2k F r)/r 3 , that originates from electron 
backscattering with the change of momentum equal to 
2k F , twice the Fermi momentum^. Correspondingly, in 
a conventional two-dimensional electron gas 3 - the ampli- 
tude of these Friedel oscillations decays as only the sec- 
ond power of the distance oc r~ 2 . In systems with the 
Berry phase and non-trivial chiral spectrum Friedel oscil- 
lations could still decay faster, oc r -3 , if the states with 
momenta k and — k are orthogonal and backscattering 
is suppressed^. In particular, this happens in case of a 
doped graphene^, when the Fermi level is shifted away 
from the Dirac points, k F ^ 0, with pscudospin-diagonal 
impurity potential. For a short-range impurity with the 
potential V(r) = US(r) the resulting induced electron 
density is&£ (we use units with H = 1) 



n(r) 



U 



27T 2 1T' 



■ cos(2fc F r), k F r ^> 1, 



(1) 



here v is the velocity of Dirac electrons. Within the lin- 
ear response approximation the derivation of this result 
is straightforward with the help of the density-density 
correlation function II w (r — r') taken in the static tu = 
limit: n = IIq (8 V. In case of a short-range potential the 
latter formula gives simply n(r) = UHo(r), and eventu- 
ally leads to Eq. flTJ. 

The calculation outlined above can be extended easily 
to the case of intrinsic graphene, k F = 0, or, cquivalcntly, 
to short distances, r <C k F x . Interestingly, the functional 
form H]) is recovered again, as can be verified within the 
same linear response approach; the only change being in 
the numerical factor 7r/4: 



n(r) 



V 
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k F r < 1. 



(2) 



The non-integrablc singularity at r — > makes the to- 
tal induced electron density J n(r)d 2 r diverge as a power 



law. This indicates the failure of the first Born approxi- 
mation when U/vr ~ 1. Fortunately, the delta- function 
potential allows for an exact non-perturbative solution 
via the T-matrix in terms of the electron Green's func- 
tion, T(E) = U/[1- W£ v G(E,p)}, which yields, 



T(E) 
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The ultraviolet divergence is cut-off at short distances of 
the order of the lattice spacing a. 

In the "unitary limit" of strong interaction, U — > oo, 
the resonant form of the T-matrix © clearly should lead 
to a significant modification of the dependence of the in- 
duced density on distance. It is not difficult to obtain a 
rough estimate for the effect: since typical energies are 
\E\ ~ v/r, we immediately conclude that in the logarith- 
mic approximation, n(r) ~ l/[r 2 In (r/a)]. Surprisingly, 
as we demonstrate in the present paper this crude guess 
is incorrect and the induced density in the unitary limit 
is in fact suppressed according to (U > 0) 



n(r) 



7rr 2 In (U /va) 



(4) 



The behavior described by Eq. (j4]) cannot be properly 
accounted for if the imaginary part of T{E) is neglected. 
It turns out that despite the imaginary part in the de- 
nominator of Eq. ([3]) being small compared with the log- 
arithmically large real part, the two parts lead to the 
contributions that largely cancel each other. 

The suppression of the local density (@| at the unitary 
limit does not mean that the total density vanishes too. 
Quite to the contrary, as the range of distances, r <C U/v, 
where Eq. (U) is applicable, increases with increasing U 
the total induced density tends to a limit / n(r)d 2 r —> 4, 
i.e. one electron per spin/sublattice. 

The resonant 1 / _E-behavior of the scattering amplitude 
in graphene has significant implications. In particular, 
the unitary limit leads to singular corrections to the low- 
energy density of states^— . The effective Casimir-likc 
coupling between two point impurities (adatoms) is pre- 
dicted to become long-range in the unitary limitii with 
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the potential energy oc [rln(r/a)] _1 . However, achieving 
the unitary limit by the strength of the potential alone 
might be difficult. To illustrate this for atomic impu- 
rities let us consider the simplest microscopic model of 
a point impurity that yields the expression ([3]), namely, 
the tight-binding approximation for the honeycomb lat- 
tice where a single carbon atom is substituted with an 
impurity ato m 12 ' 13 . In that case the effective strength of 
the delta- function is simply U ~ (V)a 2 , where a = 1.4 A 
is the interatomic distance and (V^) is the expectation 
value of the impurity's potential energy (calculated with 
the help of the unperturbed orbital). We therefore esti- 
mate that for the unitary limit to be effectively achieved 
at distances r the strength of the impurity should exceed 



(Vi) » t 



r/a 



In (r/a) 



(5) 



with t = v/a 3 eV being the hopping energy. From 
the expression ((SJ) it follows that already at distances 
r ~ 1 — 2 nm the impurity strength (Vi) needs to be of 
the order of tens of eV, which is clearly impractical. 

This restriction could be lifted if an impurity has a lo- 
calized level with low energy £o "C t leading to a resonant 
enhancement of scattering^. The limit of strong impu- 
rities could also be realized in vacancies^. Yet another 
possibility is a finite-size scattere r 16 ' 17 . Indeed, consider 
an impurity potential of finite radius p. Treating po- 
tential energy pcrturbatively is justified as long as the 
characteristic potential energy Vo is small compared with 
the typical kinetic energy, v/p. The latter estimate uses 
that the typical electron momenta inside the potential in 
intrinsic graphene are determined by the width of the po- 
tential, k ~ 1/p. The Born approximation, Eqs. ([TJ-CO), 
then holds true for small effective dimensionless coupling 
constants, g = Vop/v <C 1. The non-pcrturbativc regime 
is reached when g ~ 1, which is a much less stringent (by 
a factor r/a) condition than the above discussed condi- 
tion for a delta-function impurity. This is also evident 
from the correspondence irVgp 2 — > U expected to exist 
between the potentials at p — > 0. 

In the present paper we address the response of Dirac 
electrons in graphene to strong impurities. Since the case 
of a finite-size scatterer reveals the richest behavior and 
is free from low-distance singularities we are going to ad- 
dress it first. Still, the results obtained will be applicable 
to other strong impurities with proper modifications. In 
Section II we begin with analyzing the square potential 
problem, V^r) = V)Q(p — r), and finding the induced 
density in case of intrinsic graphene. The unitary limit is 
resonantly achieved whenever the coupling constant co- 
incides with a zero of a Bcssel function, Jo(g) = 0, where 
the induced density changes its dependence on the dis- 
tance to oc 1/r 2 . Physically, resonances are related to the 
number of wavelengths induced by the potential Vo that 
fit inside it. Such a possibility is absent for point-like 
scatterers where the strength of the potential is the only 
variable parameter that has to be sent to infinity in or- 
der to achieve the unitary limit. In Section III we address 



the effects that occur due to finite Fermi momentum in 
doped graphene, including Friedel oscillations. 



II. IMPURITY-INDUCED ELECTRON 
DENSITY IN INTRINSIC GRAPHENE 

Solutions of the Dirac equation in a centrally symmet- 
ric potential V(r) have been considered extensively be- 
13lfl6lll8l . but for the sake of convenience 



fore since Refs. 
we present them here again. The Dirac equation 



d x - id,, 




y )yj = [E- V(r)]iP, (6) 



d x + id y 



determines the two-component wavefunction is a linear 
combination of partial waves with the angular momen- 
tum quantum number m, = ±1/2, ±3/2, .. 
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Here we introduced the polar coordinates: x + iy = re 1 ^ . 

We first focus on the repulsive potential Vo > 0. As 
the integral over all filled states (negative energies E) has 
to be taken eventually, we parameterize E — —vk with 
positive k. Then the Schrodingcr equation (J6j> gives two 
coupled equations for r < p, 

m - 1/2 .A., Vb\ 

d r (pi tpi = -i [k H ) ip 2 , 

r \ v J 

m + 1/2 ■f h ,Vo\ , a , 

d r ip2 H v?2 — * k H (pi. (8) 

r \ v J 

Similarly, for r > p the wave functions obey the same 
equations as Eqs. ((HJ) but without Vq. Both inside and 
outside of the potential the wave functions satisfying the 
condition of regularity at r = are given in terms of 
Bessel functions, 



¥>2 



= A 



J m -i/2 (kr + gr/p) 
-iJm+i/2 ( kr + 9 r /p) 



(9) 



B 



C 



sgn(fc)y m _ 1/2 (|fc|r) 
(\k\r) 



-iY, 



m+l/2 



for r < p, and 

sgn(fc)J m _i/ 2 (|fe|r) 
-iJ m +i/2(\k\r) 

(10) 

for r > p. The coefficients A, B, and C arc found from 
the conditions of continuity at r = p. The normalized 19 
wave functions for m = 1/2 in the outer region, r > p, 
are then found to be 



tp! = sgn(fc) 



<P2 = -I 




/3J (|fc|r)+ 7 y (|fc|r) 



PM\k\r) + 7 n(|fc|r) , (11) 
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where the following coefficients are defined (k > 0), 
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9)Yo(kp) 
g) Ji(kp) 
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Ji (kp- 



g)Yi(kp), 

g)Jo(kp). (12) 



The electron density n = J dk(\tpi | 2 + |<^2| 2 ) takes into 
account the contributions from the two leading channels 
with m = ±1/2 as well as spin/ valley degeneracy. Sub- 
tracting the equilibrium (g = 0) density we find, 

n(r) = J ^f{ 2 ^ Mkr)Ukr) 



jL^ J fj2 + 

i=0,l J 

+ 7 2 [Y 2 (kr)~J 2 (kr)]}, 



(13) 



The latter equation can be written more compactly if 
we use the Hankel function, H^(z) = J{z) + iY(z), to 
obtain 



n(r) = — Im 

7T 



E 

i=0,l 



dk kj 



H?\kr) 



(14) 



We are now going to analyze this expression for different 
values of the impurity strength g. The oscillating behav- 
ior of the Hankel functions ensures that the integral in 
Eq. (|T4"|) converges at k ~ 1/r. Since we are interested 
in the behavior at long distances r » p we can write for 
generic values of g (though still 3> p/r), 



i ~ -Ji(g)- 



(15) 



The latter approximation is valid unless g is very close to 
a zero of the Bessel function Ji (g) (see below) . Similarly, 
by keeping the singular terms in the small-argument ex- 
pansion of the Bessel functions we can write for the sec- 
ond coefficient, 



P « Jx(g)Y (kp) - J (g)Y 1 (kp) 



J ± (g) ln(fcp) + J Q (g) 



kp 



(16) 



Substituting approximate Eqs. (|15H16[) into the exact for- 
mula (|14[) we arrive at the induced density, 



n(r) 
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pzj J\{g) pz 



7T 
I — 

2 



(17) 



Here we introduced z — kr. The integral is most easily 
calculated by deforming the integration contour so that 
it follows the positive half of the imaginary axis. The 
Hankel functions have a branching point at z = and 
a cut that extends along the negative real axis. With 
such a choice of the branch cut the logarithm in the de- 
nominator takes the value ln(r/yp) — iit/2 on the pos- 
itive part of the imaginary axis (z ~ iy). The Hankel 



functions become Macdonald functions via 
75 (— [Ki(y)] 2 . We therefore arrive at 



H\ l \iy) 



n(r) 
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(18) 



Due to the exponentially suppressed MacDonald func- 
tions, this integral is dominated by y < 1. Then, un- 
less the combination Jo(g) / J\{g) is small, the logarithm 
in the denominator can be safely neglected. Utilizing 
the numerical value of the integral, J °° dyy 2 [KQ{y) — 
K Kv)} = -^/IG, leads to the following expression, 



n(r) 



■Mg) p 

Mg) 4r 3 ' 



(19) 



This equation determines a nontrivial g-dependence, rep- 
resented in Fig. [TJ For weak impurities, j<1, Eq. (JT5J) 
reproduces the first Born approximation for the delta- 
function ([2]), if one takes into account the obvious corre- 
spondence for the strength of the delta-function, irvgp — > 
U. For strong impurities g > 1 the induced density is not 
positively defined and can become negative, which hap- 
pens for example when 2.4 < g < 3.8, sec Fig. [1] 

When the impurity strength obeys the equation 
Ji(g) = 0, the induced density is strongly suppressed, 
e.g. for g = 3.8. In such a case it is no longer possible to 
use the approximation (|15p. From Eq. (|12[) we now get 
7 « fcp(Jo/2 — J{). The presence of an extra power of k 
ensures much faster decay of the electron density 



1 \ P 
n(r) oc — . 

r 4 



(20) 



So far we have discussed the contributions from the 
lowest order s-wave scattering, m = ±1/2. A few words 
are now in order about higher in channels. The equations 
(|!Hll2p are fully applicable there as long as the order of 
the Bessel functions is adjusted: — > m — 1/2, 1 — > 
m+ 1/2 (for positive to). Instead of Eqs. f|15til6[l we now 
have 7 oc (kp) 771 ' 1 / 2 , /3 cx \/{kp) m+1 / 2 . From Eq. (13} 
we observe that each extra order of to brings additional 
factor (kp) 2 . Upon taking the /c-intcgral we conclude that 
already the p-wave scattering contribution is suppressed 
by a small factor (p/r) 2 and could be safely neglected. 
It is also known that for \m\ > 3/2 zero energy bound 
states exist 2 ^ when J\ m \-i/2(Vop) = 0. However, as the 
electron density for these states decays oc r^ 2 !" 1 ! -1 their 
contribution to n(r) is negligible. 



A. The case of resonant scattering 

With increasing g the system passes through a set of 
resonances^ determined by the condition Jo(g c ) = 0, 
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FIG. 1: (Color online) Dependence of the induced den- 
sity (rescaled by r 2 ) on the dimensionless impurity strength, 
g — Vop/v, plotted from Eq. (|18l) for r/p = 10. The approx- 
imation which is valid away from the resonance, Eq. (|19p . is 
shown for comparison with the green dashed line. The most 
notable feature is the reversal of the sign of the induced den- 
sity from positive (electrons) to negative (holes) as g passes 
through the point where Jo(g) = 0. Near the "anti-resonance" 
point g cs 3.8, determined by Ji(g) = the impurity becomes 
"invisible" as the induced density is strongly suppressed by a 
small factor p/r. 



the first of which occurs at g c = 2.4, Fig. [TJ The phys- 
ical origin of these resonances is quite clear and could 
be identified with the number of potential-induced Fermi 
wavelengths v/Vq that fit inside the radius of the well p. 
In the vicinity of the resonance Jo(g) in Eq. (fl8|) could 
be expanded as 



Mg) ~ -Ji(9c)Sg, 



9 = 9c 



Sg. 



(21) 



Contrary to the calculations of the preceding section, the 
logarithm in the denominator of Eq. (fT5|) can not be 
neglected anymore. In fact the main contribution to the 
integral comes from small values of y. The numerator 
mostly comes from K^(y) m y~ 2 , yielding, 



n(r) 
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dy 
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In agreement with the assumption just made the ar- 
guments relevant in the integral (|22[) are small, y ~ 
Sgr/(p\ln5g\) -C 1. With the logarithmic accuracy (ne- 
glecting double logarithm) we obtain, 



2 sgn(cJg) . p (r 

n(r) = — - . . , (5 5 <-ln - 
7rH In jogf| r \p 



(23) 



FIG. 2: (Color online) Sensitivity of the resonant behavior to 
small levels of doping, UfV = 0.05, and r/p — 10. Induced 
density n(r)r 2 is plotted near g = g c + 5g. The black dashed 
line represents the behavior of induced density in undoped 
graphene. In case of p-doping the resonance is shifted by 
Ag c — k F p\n(k F p) < 0. Similarly, in case of an ?i-doped 
system the changes are reversed (for the discussion of the n- 
doped graphene see Section III.C). 



with the ultimate logarithmic suppression in the immedi- 
ate vicinity of it. The dependence of the density on dis- 
tance changes from r -3 away from the resonance to r~ 2 
near it. The maximum of the density occurs at impurity 
strength \Sg\ ~ (p/r) In (r/p) where the two expressions 
([2"31 and (TK))) match, as could be expected. The behav- 
ior of the induced density in the vicinity of resonance is 
presented in Fig. [2j 

It is now interesting to compute the total induced elec- 
tron density. Such a calculation can be most rigorously 
carried out by integrating the general expression (|18p : 
it is convenient to first integrate over r (from p to oo) 
and then over y. The same results, however, could much 
faster be found straight from Eqs. (Tr9|) and (j2"5")) . Far 
from the resonance, g <C 1, integrating Eq. (| 19[) we ob- 
tain, Ntot = \g- As expected, in this perturbative regime 
the total disturbance of the system is small. Near the 
resonance, Sg -C 1, upon integrating Eq. (|23[) from p to 
p/\Sg\ we arrive at 



Ntot = -4 Bgn(5ff). 



(24) 



Thus, the total number of particles induced by the reso- 
nant impurity is equal to one per Dirac cone (consider- 
ing spin and valley degeneracy) and changes sign right at 
the resonance. This conclusion is in agreement with the 
Fricdcl sum inleM. 



B. Delta- function potential 



This relation describes the replacement of the growth of 
the density when a resonance is approached, sec Eq. (|19l) . 



We now discuss a connection between the finite size 
well and the point impurity V(r) = U6(r). As already 
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evident from Eq. (|19[) there is no rigorous limit of pro- 
gressively narrower and higher potential, g = Vp/v cx 
1/p oo, as p tends to zero. Alternatively, this could 
also be seen from the fact that the delta-function poten- 
tial does not allow a dimensionless coupling parameter 
independent of the distance r of the kind that g is for 
a finite-size well. Yet, there is a connection between the 
two models. 

First we note that the combination of f3 and 7 in the 
integrand of Eq. (| 14|) is proportional to the T-matrix. 
Indeed, by taking the asymptotics of the wave function 
(|TT|) the scattering phase shift could be expressed as, 
tan 5 = —7//?. The relation between the phase S and 
scattering amplitude in graphene is very similar to the 
standard quantum mechanical expression^ and was writ- 
ten in Ref. 0. In our notations, 



T(k) 



4w 



1 



k /3 + i-f 



(25) 



Comparing now this general expression with Eq. ([3]) we 
can identify the phase shift for a point impurity, 



cotS =-?■ = — -~ln(— 
7 kU ir \ka 



(26) 



From Eqs. (|15[) and (| 16[) we can now extract the corre- 
spondence between the two problems, 



U 



2irvp 



■h(g) 
Ma)' 



(27) 



The combination in the right-hand side (without v) could 
be identified with the scattering length of the circular 
potential welU£. 

We observe that the correspondence already encoun- 
tered away from the resonances, cf. Eqs. @ and (THJ1) . 
holds in fact everywhere. The finite-size resonances dis- 
cussed in the previous section thus represent realizations 
of the unitary limit. Replacing now Sg — > —2nvp/U in 
Eq. (|23p we arrive at Eq. (U) up to a replacement of p 
with a under the logarithm; this distinction is beyond 
logarithmic accuracy anyway. 

It should be emphasized that the imaginary part in 
the scattering matrix Eq. ([3]) cannot be neglected. This 
can be seen following the transformation from Eq. p7|) 
to Eq. (fTS)) . where it is important that the imaginary 
part is cancelled by the phase coming from the logarithm 
upon the rotation to the imaginary axis. Without in/2 
in the denominator of Eq. (fTT)) . the density would have 
been finite in the unitary limit instead of the suppression 
described by Eq. (01 . 

Note that the correspondence (|27|) is different from 
what one would expect to be the true delta-functional 
limit: U -<->• nvgp with p — > 0, g — > 00. Interestingly, 
the two models coincide with each other only in the case 
of a shallow well, j « 1. This situation is unique for 
the linear spectrum of graphene. Indeed, in the case of 
a conventional parabolic spectrum, e = p 2 /2m* the di- 
mensionless phase shift 5 = 5(Vop 2 m*) could only be 



a function of the combination Vop 2 that becomes U/ir 
in the limit of a deep and narrow well. This could be 
qualitatively understood by recalling that the parameter 
Um* gives the number of bound states in a well of depth 
Vo = U/irp 2 , which does not depend on Vq or p sepa- 
rately. In graphene, to the contrary, the phase shift is a 
function of g = Vop/v. 

A closely related to the delta-function potential is 
the case of point-like Anderson impurity with a low- 
energy state £0 hybridized with graphene ban electrons 
via Hi n t = u(dJip(0) + c.c), the problem considered in 
Ref. [HI- When £q ps the T-matrix has the resonant 
form ^ with the strength of the potential replaced by 
U —t u 2 /eo- Similar substitution in Eq. (|4]) yields the 
distribution of the induced density. Note that the density 
is positive as long as £0 > (electrons are expelled from 
the vicinity of the impurity) and negative when £0 < 0. 
Hydrogen adatoms are known to have resonances very 
close to the Dirac point^i, £ ~ 30 meV. 



III. EXTRINSIC GRAPHENE, k F 

Let us now consider the case of gated or doped 
graphene with a nonzero Fermi-momentum. We will an- 
alyze two cases separately. When the distance to the 
impurity exceeds the Fermi wavelength, k F r 3> 1, the 
conventional Friedel oscillations develop whose specific 
behavior depends strongly on the impurity strength. The 
opposite case of a weakly gated graphene, k F r <C 1, is 
much more spectacular. Because of the resonant behav- 
ior described by Eqs. (fT5|) and (JD) the sensitivity of the 
induced density to small k F turns out to be very strong. 



A. Weakly p-doped or gated graphene, fepr <C 1 

The case of a p-doped graphene is particularly straight- 
forward as it is sufficient to simply replace the lower limit 
in the integral in Eq. (jTTJ) with k F r. The integral is 

then written as — J^ FT ', with the first one yielding 
the same expression as before, Eq. (|2|). The second inte- 
gral brings the correction, 



k F r 



An(r) 



2 

nr 



zdz 



z In ( — 

pz 



(28) 



From the form of this integral it is clear that the correc- 
tion is highly asymmetric: it quickly decays with increas- 
ing 5g when the latter is positive. For negative Sg < 
the correction (|28|) is much stronger as the integrand has 
a pole near the real axis. The pole is significant as long 
as 2k F p\n(k F p) < Sg < 0. Outside of this region the 
correction again decreases quickly. This is illustrated in 
Fig-Hi the plot of the total density n to t(r) = n(r)+An(r) 
is effectively "shifted" to the left of g c . The first terms 
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in the Taylor expansion in powers of Sg could be easily 
extracted: 



n to t{r) 



irr 2 In {k F p) 



1 



(tt - 2)8g 
k F p\n 2 (k F p) 



(29) 



The shift of the position of the resonance due to dop- 
ing, Ag c , with great precision is given by the condition, 
j3(k F ) = 0, yielding 



Ag c = — fc F pln 



1 

k F p 



(30) 



see also discussion in the next Section. Another conse- 
quence of Eq. (|29]l is that the singularity in the derivative 
of the undoped resonant density is regularized by finite 

Let us emphasize the signature feature of the resonant 
behavior depicted in Fig. [3J namely the reversal of the 
sign of the density via small doping occurring in the range 
Ag c < Sg < 0. 



1. Delta- function potential 

While the effects of weak doping lead to significant 
qualitative modification to the dependence of induced 
density on the distance for finite-size impurities, the cor- 
responding changes for point defects are even more dras- 
tic. Note that the induced density in the latter case is 
always of the same sign as U in case of intrinsic graphene. 
In other words the delta-function potential roughly maps 
on the interval < g < 2.4, below the first resonance of 
the finite-size well, which corresponds to the U — >• oo 
unitary limit, cf. Eq. (|27j) . Doping with small amounts 
of holes (for U > 0) or electrons (for U < 0), on the 
other hand, pushes the resonance to a finite impurity 
strength Uq, so that for U > Uq the induced particle 
density changes sign, much like in the case of a finite-size 
well. This strong modification of the induced density is 
captured in Fig. [3l 

Although we explained below Eq. (|27p that obtaining 
expressions for U S(r) is rather simple, let us present here 
for reference the total induced density for the latter case, 
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(31) 



where by u = U/vr we denoted the dimcnsionlcss param- 
eter of the problem. Fig. [3] illustrates the dependence of 
the induced density on u for different values of doping 
levels kpr. The density changes sign at the point where 
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(32) 
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FIG. 3: (Color online) The emergence of a resonance for a 
point-like impurity described by the delta-function potential 
V(r) — US(r) for weak p-doping levels (for U > 0) versus 
the dimensionless impurity strength u = U /vr, plotted from 
Eq. Q3ip , for r/a — 10. The green dashed line illustrates 
the undoped case: the logarithmic decay Q is very weak on 
the scale shown. Finite k F makes the density change sign at 
Uo given by Eq. (|32p . The same resonance is revealed via 
n-doping for an attractive potential of the impurity, U < 0. 



At U = Uo the scattering phase shift at the Fermi level 
5 F , determined by Eq. (|25p undergoes a jump from 7r/2 
to -tt/2. 



B. Strongly p-doped case, fcpr ^> 1 

With further increase in the doping level the vicinity of 
the Fermi surface begins to dominate. Eq. (fT4")l with the 
lower limit replaced by kp r is a convenient starting point. 
Using the asymptotic expression EtoiP* Wl 2 = 

— exp (2iz) and noting again that the ratio j/f3 — 

— tan d is related to the phase shift 5 we obtain, 



n(r) 



-Im 



dz 



v6e 2%z+%s . 



(33) 



k F r 



Because of the strongly oscillating behavior of the inte- 
grand only the vicinity of the lower limit contributes to 
the density. After simple integration we arrive at, 



n(r) = 2 sm5 F 



cos 



(2k F r + S F ) 



ir 2 k F r 3 



(34) 



with the scattering phase shift taken at the Fermi sur- 
face, S F = -tan _1 [7(/c F )//3(fc F )]. The formula (JM|) de- 
scribes Friedel oscillations both in the perturbative limit 
where 5 F <C 1 and Eq. |T]) is recovered, and the resonant 
(unitary) regime, when 8 F is close to ±7r/2 where the 
amplitude of the oscillations no longer depends on the 



7 



impurity strength. Note that in the vicinity of the res- 
onance the density reverses its sign, similar to the 
intrinsic case. 



a resonance, the phase shift jumps by 7r and the sign of 
the induced density reverses, much like in Fig. [T] 
IV. SUMMARY AND CONCLUSIONS 



C. ?i-doped graphene 

Now we discuss n-doped graphene with ground state 
filled up to a positive Fermi energy, E F > 0. According to 
our notation, k = —E/v, in the n-doped case the range of 
positive energies with k F < k < appears in addition to 
the filled lower Dirac cone. For negative k the eigenfunc- 
tions Eq. ([fT]) are still applicable with coefficients now 
assuming the form 

P = J X (kp + g) Y (\k\p) + J (kp + g) Yidklp), 

7 = -J (kp + g) Ji(\k\p) - J t {kp + g) J (\k\p), (35) 

implying that Eq. (|15j) holds for 7 with both signs of 
k, whereas Eq. (TO)]) captures correct (3 for all values of 
k if one changes k to |fc| in the argument of logarithm. 
Taking this into account, the induced density is given by 
Eq. (fT3|) or (fl4|) . provided that the lower limit of the 
integral is extended to k F and k is changed to |fc| in the 
rest of the integrand. We write it as 



i(r) = -Im V 

7T ^ * 




H?\\k\r) 



(36) 

It is straightforward to check that the first integral here is 
zero. Indeed, as a combination of integrals over (—00, 0) 
and (0, 00), it amounts to the sum of two integrals, one 
of which is the same as in Eq. (|17|) . while the other one is 
given by Eq. (fT7|) with the opposite sign of Jo(g) / Ji(<?) 
in the denominator. According to Eq. (fT5j) such a com- 
bination is zero. 

We conclude that the induced density is governed by 
the second integral of Eq. ([36]) . leading to relation (j34]l 
with the opposite sign. Therefore, the general relation for 
induced density (Friedel oscillations) in strongly doped 
graphene is 



n(r) 



2vsind F cos(2|fc F |r + S F 

eZ ^ 



(37) 



The phase shift at the Fermi surface is determined by 
(cf. Eq. [261) 



cot S f 



2v Mg) 



■wE F pJi(g) 7r \\Ep\p) 



(38) 



The expression (|3T|) generalizes Eq. (TT]) and has the latter 
as the limiting case at weak couplings. The long-behavior 
is always oc r -3 , even at resonance, as opposed to the 
intrinsic graphene, Eq. (U). Still, as one passes through 



In this paper we considered the two non-pcrturbative 
models for impurities in graphene: a substitution atom 
described by V(r) = U6(r), and a finite-size impurity 
(molecule or nanoparticle) with V(r) = Vb©(p — r). The 
first model cannot be derived from the second one via 
a standard limiting procedure in the effective low-energy 
Dirac fermion description, although there is a simple cor- 
respondence between the two. As an illustration of the 
scattering problem we calculated the impurity-induced 
density both for intrinsic and extrinsic (doped) graphene. 

The case of a finite-size potential reveals a set of res- 
onances that correspond to the zeros of the Bessel func- 
tion, Jo (g) = 0, for the dimcnsionless impurity strength, 
g = Vqp/v, see Fig. [1] Near a resonance (which is an 
effective realization of the unitary limit) the r~ 3 de- 
cay of the induced density is replaced with a slower 
r~ -dependence. The resonances are interspaced with 
"anti-resonances" that occur when J\(g) = 0. At the 
latter points the induced density becomes strongly sup- 
pressed with the long-distance behavior oc r -4 . Low dop- 
ings, kpr <C 1, modify the vicinity of resonances lead- 
ing to the shift of the resonant coupling g c and removal 
of the logarithmic singularity, Fig. [2] The rest of the 
n(g, r)-dependence is unaffected by low k F , including po- 
sitions of the anti-resonances. At stronger doping levels, 
kpr ^> 1, the induced density follows the usual Friedel 
cos (2kpr) /?* 3 -dcpcndcnce known from the first Born ap- 
proximation, with the resonant behavior entering via the 
phase shifts. In particular, near the resonances the in- 
duced density reverses sign. The Friedel oscillations re- 
main strongly suppressed near anti-resonances. 

The point impurity model V = US(r) does not have 
anti-resonances. In the intrinsic case the only true res- 
onance occurs at U — 00 though the maximum of the 
induced density is found near U ~ lOvr. With the dop- 
ing the resonance is shifted to the finite values Uq given 
by Eq. ((321), see Fig. [3 

Our findings indicate that the unitary limit of strong 
impurities could be realized with realistic potentials Vq. 
In particular, for a nanometer-size impurity p ~ I nm 
the first resonance occurs when Vq = 2Avh/p w 1.5 eV. 
Correspondingly, for an extended p ss 10 nm defect the 
potential would have to be only Vq ~ 150 meV. 
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